Friedel oscillations for interacting fermions in one dimension 
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We study Friedel oscillations in one-dimensional electron liquids for arbitrary electron-electron interaction and arbitrary 
impurity strength. Explicit results for spinless as well as spin-i electrons are given. In the case of Luttinger-liquid leads, the 
Friedel oscillations decay as far away from the impurity where g is the interaction constant. For a weak scatterer, a slower 
decay is found at small-to-intermediate distances from the impurity, with a crossover to the asymptotic decay. 

PACS numbers: 72.10.-d, 73.40. Gk 



Quasi one-dimensional (ID) interacting fermion sys- 
tems have attracted renewed attention recently, partly 
due to the technological relevance of such systems. They 
are also appealing because of the existence of exact solu- 
tions for simple models |l|^ . The striking non- Fermi liq- 
uid behavior of ID electrons is most clearly exhibited in 
transport properties, especially in the presence of impuri- 
ties or barriers Q . Usually, either the case of a very weak 
or an almost insulating barrier have been studied explic- 
itly. The crossover between these two regimes, however, 
has rarely been looked at despite its importance for sev- 
eral fundamental issues, e.g. pinning of a Wigner crystal 
I pof or charge density waves [|ll|, breakdown of charging 
effects with increasing tunnel conductance ||l^ , or trans- 
port in ID quantum wires or heterostructure channels 
p^ . In this Letter, we discuss several aspects of this 
crossover, in particular properties of the screening cloud 
around the impurity. 

The presence of an impurity in a metal is known to 
cause Friedel oscillations in the density profile due to the 
sharp Fermi surface |p^ , ^ , 

p{r) cos(2/cF?' + S)/r'^ , 

where is the Fermi vector, d the dimensionality, r the 
distance from the impurity, and 6 a phase shift. This 
result holds far away from the impurity and for a Fermi 
liquid only, and in ID the question arises how Friedel os- 
cillations are affected by the presence of strong electron 
correlations. As pointed out by Matveev et al. ||^, for 
weakly interacting electrons in ID, one can understand 
the zero- voltage anomaly and power-law temperature de- 
pendence of the nonlinear conductance |^] in terms of 
electron backscattering by Friedel oscillations. It is thus 
surprising how little attention has been devoted to the 
modification of Friedel oscillations by strong interactions. 

We treat the ID interacting electron liquid in the 
framework of standard bosonization |^ ^ . This approach 
is appropriate for low temperatures where only excita- 
tions near the Fermi surface are relevant. The electron 
creation operator for spin s = ± at position x is expressed 
in terms of the boson fields Of^{x) and (j)p{x) {fi = p, (t) 
which arise as linear combinations of spin-up and spin- 
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X exp ^i\/TT/2[(j)p{x) + s(j)^{x)] 
The boson fields obey equal-time commutation relations 

[(t>f_,{x),9^{x')] = -{i/2)Sf,^sgn{x - x') , 
and the canonical momentum for the 9^ field is 11^ = 

We are concerned with density distributions in the 
presence of impurities or barriers. The bosonized form 
of the density operator is |^ 



p(x) = yy2/Trda:0p{x) H cos[2/cFa; + v2^6'p(a;)] 

TT 

X cos[A/27r6'CT(a;)] + const. cos[4fcFa; + V8Tr0p{x)] , (2) 

where the background charge po — 2k-p/'rT has been 
omitted. The three terms in Eq.(^ are (1) the long- 
wavelength contribution, (2) the 2k-p charge density wave 
part, and (3) the 4fcF Wigner component |^J^. The 
Wigner component is not present in the spinless case, 
since two right-movers have to be flipped into left-movers 
simultaneously for this term to arise. 

Assuming to be away from lattice or spin density wave 
instabilities, and neglecting electron-electron backscat- 
tering for the moment, the clean system is described by 
Hn = H„ 



Hp with 



Hp=.^Jdx [Ulix] + idx9pix)r] 
dxdx' U{x — x') dxdp{x) dx 



(3) 
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where is the Fermi velocity and Ti ^ 1. The spin part 
H„ is identical to the charge part Hp with no interaction 
potential U and the p fields replaced by the a fields. 
Here, U{x) is a (screened) Coulomb interaction, and we 
will explicitly study a short-ranged potential (Luttinger 
liquid) [pi and a 1/r long-ranged potential Isir 
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Let us now consider a scattering potential. Assuming 
an essentially pointlike scatterer at a; = 0, one finds 



i/top = ycos[x/2^0p(O)] cos[V2^(?„(0)] . 



(4) 



Spin and charge parts are now coupled through this term 
1^. Actual computations using the bosonized model 
Hhos = ^^0 + ^^imp necessitate introduction of a cutoff 
parameter lu^ | p6[ . Increasing V/lOc from zero to infin- 
ity corresponds to tuning the barrier from transmittance 
one down to zero. Near zero transmittance, the weak- 
link model used in Ref. |^ is reproduced by the instanton 
treatment of i?bos- One may also show by direct com- 
parison with the exactly solvable Fano-Anderson model 
that -ffbos can reproduce the full crossover quantitatively, 
thus validating Eq.(||) for arbitrary V/cOc- Although this 
comparison can be carried out only in the absence of 
Coulomb interactions, it implies that a complete descrip- 
tion of the crossover is indeed possible using i/bos- 

Friedel oscillations can be extracted from the generat- 
ing functional 



exp 



(5) 



where the average is taken over H^os- Since the impu- 
rity influences 9^ only at a; = 0, we constrain 6^{x — 0) 
to be equal to new fields, say, = ^/2tt9^{0). Repre- 
senting these constraints by a Fourier functional integral, 
one can integrate out all 0^ {x) modes due to their Gaus- 
sian nature. In contrast to previous treatments of this 
problem, we keep explicit information about the electron 
liquid away from the barrier. The remaining auxiliary 
field integrations coming from the above constraints are 
Gaussian and hence also performed easily. In the end, 
one is left with the nontrivial average over the fields 
alone, which are coupled to each other through -ffimp- 
Collecting together all terms, we obtain 



Z = B{x,X^,V) Yl W^\a 



(6) 



The functions W^(x) are independent of the barrier 
height since they do not participate in the q average, 

Wf^ix) = exp ( - 2^ ) , (7) 
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where (3 — l/k-oT and 
F^P'){x)^vv I dk- 



cos(fcx) 



(8) 



Here, LUn = 2Tm/P are the Matsubara frequencies and Uk 
is the Fourier transformed electron-electron interaction 
|l8|. A Luttinger liquid is governed by the interaction 



constants = 1 and gp = g < I (|J^. In that case, 
Eq.(@) becomes simply 



^;l^)(x) = f^exp 



\g,i(^nX\ 
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Finally, the quantity B in Eq.(|6|) is an average in q space; 
all dependency on impurity properties is contained in this 
factor. With Matsubara components q^.n, we find 



iX,, 



Fi^\x) 



J / q 



The q bracket stands for an average taking the action 
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V / dTC0s[5^(T)] cos[(7p(t)] 



The first term is the standard influence functional 
From these equations [or generalizations with additional 
fields at some other position x'] , one may obtain all de- 
sired information about density profiles and correlations 
in presence of an arbitrarily high barrier. The results of 
Refs. [|[|| for the clean system {V — 0) are easily re- 
covered. Similar expressions incorporating the cj) fields 
reproduce the results of Refs. |^,^. 

Let us now consider the expectation value {p{x)). From 
Eq.(^, one finds (0^(a;)) = 0. The long- wavelength part 
in p [first term in Eq. (|^)] does not feel the impurity since 
^imp does not contain forward scattering terms. How- 
ever, Friedel oscillations follow for the 2A:f and (in the 
spin-^ case) for the 4/cf component in Eq.(^. 

For a spinless Luttinger liquid, we obtain the Friedel 
oscillation 

{pix))/po = -P{\xlg,V)Wi\xlg) cos(2fcFx) , (9) 

where pa = k-p/n. Evaluation of B gives the pinning 
function 



P{x,g,V) 



^ OO 



,-gx|w„|/t>F 



(10) 



This function determines the amplitude of the Friedel 
oscillation and hence the ability of the scatterer to pin 
charge density waves. 

In the following, we discuss the ground-state properties 
of the Friedel oscillation (|^) in some detail. From Eq. (^ , 
we find 



W{x,g) = {l + x/ay 



(11) 



where a = VF/2gujc is a microscopic lengthscale, say, a 
lattice spacing. The properties of the pinning function P 
can be studied using either numerically exact quantum 
Monte Carlo (QMC) simulations or simple approxima- 
tions. For transmittance one (V^ = 0), the "charge" q is 
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free and P = 0. For zero transmittance (V oo), the 
potential Vcosq locks q at odd multiples of tt, and P 
takes its maximal value, P = 1, for all x. 

To estimate P for arbitrary V/uJc, we first discuss 
a simple variational procedure based on a quadratic 
trial Hamiltonian (self-consistent harmonic approxima- 
tion, SCHA) |19[2^. Replacing the cosine term by a 
Gaussian with frequency Q, Feynman's variational prin- 
ciple leads to the self-consistency equation 



n = v 1 



(12) 



Within the SCHA, Eq.(|l^ is a Gaussian average, and 
one finds 

P{x,g,V) = exp[-ge^^+''y^° Ei{{x + a)/xo)] (13) 

with the exponential integral Ei(jj) and the crossover 
scale 



Xo 



27^1} 



(14) 



In the strong-scattering limit, ttV/lJc '3> 1, Eq.(p^) 
yields fl — V. In this limit, only small fiuctuations 
around the minima of the cosine potential are possible, 
with interwell tunneling being forbidden by an exponen- 
tially small WKB factor. Since xq is even smaller than 
a, see Eq.(^^, the term "crossover" is not meaningful in 
this limit. Using asymptotic properties of Ei{y), Eq.(|lJ 
becomes for x ^ a 



(15) 



in accordance with our QMC results and a recent study of 
Friedel oscillations by open boundary bosonization pll |. 

In the weak-scattering limit, ttV/uJc <C 1, the pinning 
function exhibits more structure. From Eq.(|l2|) one has 



which together with Eq.(|14|) implies that the crossover 
scale goes to infinity as ^ 0, namely xq ^ y~i/(i~9). 
For X ^ Xq, SCHA always gives P ~ 1 according to 
Eq. (|l5|) . This failure is due to the complete neglect of 
interwell tunneling in the SCHA, as can be seen by con- 
sidering the X — > oo value of the pinning function (10), 
Poo = — {cos q)q, where q is the time average value of 
the imaginary-time path q{T). Without tunneling tran- 
sitions q is an odd multiple of tt and one finds Poo = 1 
as predicted by SCHA. However, taking into account ex- 
cursions to neighboring wells, it is readily seen that in 
general Poo < 1- Despite of these shortcomings, the effec- 
tive Gaussian treatment indicates that for weak scatter- 
ers there is a crossover, with a slower decay of the Friedel 
oscillation at intermediate distances than the asymptotic 
decay. In fact, Eq.(p^) predicts P ~ a;^ for a; <C xq. 



To investigate the weak-scattering limit further, we 
have evaluated Eq.(|o|) in powers of V giving to lowest 
order 

P{x,g,V)^jg— + OiV') (16) 



with jg = {4:9-'^/n)B{l/2,g-l/2), where B{x,y) is the 
Beta function ||2^ . This perturbative result is only valid 
for g > 1/2 (otherwise jg diverges). Furthermore, since 
higher orders of the perturbation series grow faster ^ 
2;"(i-s) with n = 3, 5, . . ., the lowest-order result ( |l6| ) is 
only valid for x ^ xq where xq is found to be given by 
the SCHA crossover scale (|lj). Hence, for intermediate 
distance from the barrier, the Friedel oscillation decays 
slower than x~^ , namely like As a consequence, 

there is a nontrivial limit for the pinning function as a; — > 
oo and ^ 0. 
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FIG. 1. QMC results for the pinning function P{y) at 
g = 2/3 and two barrier heights: V/uc — 0.05 (solid curve) 
and 0.1 (dashed curve). The dimensionless space variable is 
y = uocx/ve- Note the logarithmic scales. 

The behavior of the pinning function for arbitrary 
V/lo,, and g can be computed by means of QMC sim- 
ulations. In Fig. 1, we show results for g = 2/3 and 
two (relatively small) barrier heights V . For small-to- 
intermediate x, our data display a power law P x^" 
with (52/3 = 0.24±0.03. This is in crude accordance with 
the perturbational result 5g = \ — g valid for weak in- 
teractions. On the other hand, for x ^ xq, the pinning 
function P is essentially constant, and the asymptotic 
decay of the Friedel oscillation is therefore ~ x~^ . A 
similar behavior is found aX g — 1/3 where direct per- 
turbation theory is inapplicable. Fig. 2 shows QMC re- 
sults for the pinning function. The small-to-intermediate 
X behavior is again a power law, now with exponent 
<^i/3 = 0.22 ± 0.03. This is in crude accordance with 
the SCHA prediction 5g = g which holds for very strong 
interactions. Based on our numerical data, the exponents 
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5g are independent of the barrier height while the region 
where the intermediate decay is seen shrinks rapidly as 
V grows. 
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FIG. 2. QMC results for the pinning function P{y) at 
g — 1/3 and two barrier heights: V/cJc — 0.01 (solid curve) 
and 0.05 (dashed curve). 

For spin-i electrons, we find asymptotically a slightly 
faster decay ^ a;~(i+£')/2 at zero temperature. This can 
be rationalized by noticing that the additional spin chan- 
nel has = 1 and the exponent —g in Eq.(|9|) has to be 
replaced by —{ga +5p)/2. Remarkably, for spin-i elec- 
trons, there is also a Akp Friedel oscillation component 

(p(x)) cos(4fcFa;) , 

which dominates over the 2kF contribution for strong 
enough correlations, g < 1/3. Since Akp corresponds to 
the interparticle spacing, this suggests that for g < 1/3 
signatures of Wigner crystal behavior are induced by the 
impurity. 

Wigner crystal behavior has also been found by Schulz 
1^ for the clean system with long-ranged 1/r correla- 
tions. For 1/r interactions, the Akp Friedel oscillation 
decay is extremely slow. While the spin degrees of 
freedom involve again the a;~^/^ factor suppressing the 
2A:f component, the Akp Friedel oscillations decay like 
exp(— cVlnx), i.e. slower than any power law. Effec- 
tively, one will then only observe the Akp component. 
In the spinless case, the same quasi long-ranged behavior 
appears for the 2fcF component already because the spin 
channel is absent now. 

Apparently, Friedel oscillations are always present in 
ID for arbitrary electron-electron interaction. Moreover, 
due to reduced screening in low dimensions, their decay 
is always slower than the Fermi liquid 1/x prediction. 
We wish to stress that inclusion of backscattering is not 
expected to alter these findings substantially. In the spin- 
less case, backscattering is treated as exchange event of 



forward scattering and can be absorbed by a redefinition 
of g. The asymptotic decay of the Friedel oscillation is 
then always . In the spin-i case, based on the renor- 
malization group analysis ||^ , there are at most weak log- 
arithmic corrections. 

To conclude, we have computed the Friedel oscillations 
in an interacting ID electron liquid. These results should 
show up in NMR experiments or as strong quasi-Bragg 
peaks in x-ray scattering. They are also of relevance for 
studies of quasi-one-dimensional conductors at low dop- 
ing concentrations, or the case of a magnetic impurity. 

We wish to thank P. Riseborough, M. Sassetti and U. 
Weiss for useful discussions. 
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